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ABSTRACT 

We report first multicolor polarimetric measurements {UBV bands) for the hot Jupiters HD 189733b and con- 
firm our previously reported detection of polarization in the B band (Berdyugina et al. 2008). The wavelength 
dependence of polarization indicates the dominance of Rayleigh scattering with a peak in the blue B and U 
bands of ^10~ 4 ± 10~ 5 and at least a factor of two lower signal in the V band. The Rayleigh-like wavelength 
dependence, detected also in the transmitted light during transits, implies a rapid decrease of the polarization 
signal toward longer wavelengths. Therefore, the nondetection by Wiktorowicz (2009), based on a measure- 
ment integrated within a broad passband covering the V band and partly B and R bands, is inconclusive and 
consistent with our detection in B. We discuss possible sources of the polarization and demonstrate that effects 
of incomplete cancellation of stellar limb polarization due to starspots or tidal perturbations are negligible as 
compared to scattering polarization in the planetary atmosphere. We compare the observations with a Rayleigh- 
Lambert model and determine effective radii and geometrical albedos for different wavelengths. We find a 
close similarity of the wavelength dependent geometrical albedo with that of the Neptune atmosphere, which 
is known to be strongly influenced by Rayleigh and Raman scattering. Our result establishes polarimetry as a 
reliable means for directly studying exoplanetary atmospheres. 

Subject headings: planetary systems — polarization — stars: individual (HD 189733) 



1. INTRODUCTION 

Direct detection of exoplanets with polarimetry opens 
the prospect for probing their atmospheres. A reflecting 
planet breaks the symmetry of the stellar radiation and 
conspicuously mark s its presence in polar ized light (e.g., 
iFTuri & Berdyuginall2010t lBerdyuginall201 ll) . Thanks to the 
differential nature of polarimetry, its dynamic range exceeds 
that of any other technique, but it is still a challenge to detect 
low signals. 

It is known that polarization pro perties of the light are gen- 
erally wavelength dependent (e.g., Stenflo 2005), with Thom- 
son scattering being one exception. Rayleigh scattering be- 
ing an important opacity source in upper layers of cool atmo- 
spheres results in polarization strongly increasing toward the 
blue with a A -4 law. Combined with an angular dependence 
of polarization, this dictates an optimal range of wavelengths 
and scattering angles at which polarization can be success- 
fully detected. 

There have been a few attem pts to d e tect po lari zed reflected 
light f rom an exoplanet. Hou gh etafl (120061) and lLucas et alj 
(2009, hereafter L09) observed hot Jupiters in rBoo and 
55Cnc in a red filter of 590 nm to lOOOnm (maximum at 
800 nm). Due to relatively low statistics and incomplete or- 
bital phase coverage, only standard deviations of Stokes q and 
u could be deduced : 5.1T0" 6 and 2.2-10" 6 for the two sys- 
tems, respectively. Berdyugina et al. (2008, hereafter B08) 
observed another hot Jupiter, HD 189733b, which is twice as 
close to the parent star as, e.g., r Boo b. The Stokes parame- 
ters were measured in the Johnson B band (370-550 nm, max- 
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imum at 430 nm). Despite a relatively low accuracy of indi- 
vidual measurements, a very high statistics of the data (about 
100 nightly measurements) and an even distribution over or- 
bital phases have enabled the detection of variable, phase- 
locked polarization of maximum ^2-10~ 4 and best-fit ampli- 
tudes in Stokes <? and u of (1.5 ±O3}T0^ and (1.1± 0.4>10^, 
respectively, for binned data. Wi ktorowiczl (120091 hereafter 
W09) also observed HD 189733 but in the 400-675 nm band 
(centered at 550nm). An upper limit of 7.9T0 -5 in the po- 
larization degree (5 TO -5 in Stokes q) from only six nightly 
measurements was reported. At first glance, these measure- 
ments contradict to each other. However, as they were taken 
at different wavelengths, the conclusion made from their di- 
rect comparison was incorrect. Furthermore, only one of the 
six W09 measurements (which provided the upper limit) is at 
a phase near elongation where a significant polarization sig- 
nal is expected based on our B08 ephemeris. Therefore, the 
nondetection of W09 is inconclusive and consistent with our 
detection. 

In this paper, we show that the detected amplitudes 
and upper limits across the optical spectrum indicate the 
dominance of Rayleigh scattering in the atmospheres of 
hot Jupiters. This is also favored b y transit spectroscopy 
in the optical bv iPontet alJ (120081 hereafter, P08) and 
iLecavelier des Etangs et alJ (l2008albfT and in the near infrared 
by ISing eta l. (2009). Here we report first multicolor polari- 
metric measurements in the UBV bands for HD 189733b. We 
confirm within the standard deviation the previously reported 
first detection of polarization in the B band by B08 and de- 
termine reflecting properties of the planet. We also find the 
consistency between the data of W09 and our measurements 
in the V band. 

2. OBSERVATIONS 

The new observations were carried out in 2008 April 18-24 
and August 2-9 with the TurPol UBVRI polarimeter mounted 
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in the Cassegrain focus at the 2.5m Nordic Optical Tele- 
scope (NOT). The polarimeter was designed by Piirolaj d 19731 
1988). It consists of a half- (or quarter-) wave plate as retarder, 
calcite block, diaphragms, chopper, filter wheel, dichroic mir- 
rors, and five photomultipliers as detectors for the UBVRI 
bands. The calcite block splits the light onto parallel ordi- 
nary and extraordinary beams. Thus, the sky contributes to 
both beams and its polarization (e.g., due to dust in the Earth 
atmosphere or moonlight) can be exactly compensated. This 
is a great advantage of TurPol to the polarimeters PlanetPol 
of lHoughetaH ( 120061) and its copy POLISH of W09 which 
have to monitor very accurately a variable sky polarization. 
The rapidly rotating chopper (25 Hz) alternates integration of 
the ordinary and extraordinary beams at the same detector for 
each band simultaneously, which strongly diminishes system- 
atic effects due to instrumental sensitivity, seeing and vari- 
able sky transparency. To measure linear polarization, the re- 
tarder was rotated at 22.5° intervals. Each pair of the observed 
Stokes q and u was calculated from eight exposures at dif- 
ferent orientations of the retarder, which allowed us to avoid 
systematic errors due to imperfections of the retarder. Our 
definition of the gM-system is in accordance with the common 
agreement that positive and negative q are in the north-south 
and east-west directions, respectively, while positive and neg- 
ative u are at an angle of 45° counter-clockwise from the pos- 
itive and negative q. 

To determine the instrumental polarization, we observed 26 
nearby, bright, nonvariable, and nonpeculiar stars (V=4 m -6 m , 
d < 50 pc) of spectral classes A-G, which were expected to 
have no detectable intrinsic and interstellar polarization. Note 
that W09 observed only one standard star which resulted in a 
poor absolute calibration of his data. These stars were repeat- 
edly observed at various parallactic angles to determine both 
the telescope and instrument polarization, the former rotating 
during observations due to the alt-azimuthal mount and the 
latter being constant. Some stars were observed a few times 
on different nights to ensure the stability. In addition, for 
calibration of the polarization angle zero point, we observed 
highly polarized standard stars HD 132052, HD 16 1056, and 
HD204827. Total integration time of 1-2 hours for each star 
allowed us to measure instrumental polarization with the ac- 
curacy of (1-2)- 10~ 5 . Measurements for those standard stars 
which were used for the calibration are presented in Table [T] 

The instrumental polarization deduced from these observa- 
tions and expressed in polar coordinates (the polarization de- 
gree f; ns and the angle 9 ms ) is as follows: in April 2008 
P im (U) = (23.1±2.7)-10" 5 , 6> ins (t/) = 106.5 ± 3.3, 
P ins (fi) = (13.4±1.6)-10" 5 , 6 ins (B) = 117.9 ± 3.4, 
P ins (V) = (11.8±2.5)-10" 5 , 9 ms (V) = 116.5 ± 5.9; 
and in August 2008 

P ins (U) = (24.2±1.4)-10" 5 , 9 in ,(U) = 110.2 ± 1.6, 
P ins (fi) = (13.8±1.1)-10- 5 , 6 ins (B) = 116.2 ± 2.2, 
An S (V) = (11.9±1.3)-10- 5 , ms (V) = 121.1 ± 3.1. 
We found a remarkable agreement between measurements 
carried out in April and August. This proves the stability of 
the total instrumental polarization during our observing runs 
as well as the robustness of our choice of standard stars. Note 
also that during our August run there was some presence of 
Sahara dust in the atmosphere. L09 reported that Sahara dust 
can produce sky polarization of order 10~ 5 at large zenith dis- 
tances. This is at the noise level of our measurements. Also, 
as explained above, the construction of the TurPol enables 
exact compensation of any background polarization which is 



Table 1 

Polarization of standard stars. 
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|_-o, 4J ± J 


r a on i a 
L-o, oJ ± 4 


r in 1/11 i o 
L-1U, 14J ± o 


HD1215oU 


r 1 a in I in 

[14, -ljj ± 10 


[-5, -9] ± 4 


[-9, 13] ± 6 


n \j i zouj j 




[\A ni _L s 


re ci _|_ c 
[-J, -oj rx J 


HD 157466 


[0, -7] ± 6 


[4 S -5] ± 4 


[-5, 4] ± 5 


HD174160 


[3, 0] ± 4 


[-6, -5] ± 4 


[-9, 16] ± 5 


HD 177082 


[6, -2] ± 5 


[0,-1] ±4 


[0, 4] ± 6 


HD 182807 


[-4, -2] ± 3 


[-4, -4] ± 3 


[-2, 6] ± 4 


HD 197076 


[-3, -15] ±3 


[-7, -6] ± 3 


[-2, -3] ± 3 


HD198390 


[-2, 1] ± 4 


[-5, -4] ± 4 


[-7, -6] ± 4 


HD206826 


[2, -5] ± 6 


[0, -6] ± 4 


[-8, -2] ± 5 


HD212395 


[-10,-16] ± 8 


[-5, -4] ± 3 


[-16, -1] ± 3 


HD216756 


[10, 3] ±4 


[10, -1]±5 


[9, 2] ± 4 


HD2 18261 


[1,-1] ±3 


[-1,-2] ±3 


[-10, 4] ±4 


HD225239 


[1,2] ±4 


[-2, -6] ± 3 


[1,7] ±4 



Note. — Polarization is in units of 10 5 . The upper and lower parts are for April and 
August, 2008, respectively. 

not variable within the time of alternation between the ordi- 
nary and extraordinary beams (1/25 s in our case). The very 
good agreement between the measurements of the instrumen- 
tal polarization in April and August proves that this technique 
works very well. Moreover, since we increased total integra- 
tion time in August due to increased opacity of the sky, the 
accuracy achieved in August was on average better than that 
in April. Therefore, we conclude that Sahara dust has not 
changed the polarization noise or introduced bias into our ob- 
servations. 

Our science targets were a number of nearby systems with 
hot Jupiters. Here we report results for HD 189733b only, 
while data for other systems will be analyzed in subsequent 
papers. We made 10 s exposures in the UBV passbands si- 
multaneously during 3-4 hours every clear night and obtained 
1-2 measurements per night in each band (35 in total) with in- 
dividual errors of (4-6)- 10" 5 in the U band and (2-4) 10" 5 in 
the B and V bands. The data were corrected for instrumen- 
tal polarization. The calibrated measurements for the three 
passbands are shown in Fig. Q] (upper panels). The following 
amplitudes in the UBV bands are detected (in the 10~ 5 scale): 
in Stokes q, respectively, 9.4±2.1, 9.1 ±2.4, 4.6±2.4; and in 
Stokes u 4.4±2.8, 4.3±2.4, 2.1 ±1.4. The uncertainties were 
calculated as standard deviations from the best fits. The B 
amplitude is in agreement with the previously reported polar- 
ization by B08 within the standard deviation. In Stokes q the 
largest difference is seen only in one B08 binned data point 
near the phase 0.7, while all other points agree well with the 
new data. In Stokes u only two points deviate near the phases 
0.2 and 0.65, but they are still within about one standard devi- 
ation. As the U and B amplitudes are similar, we can improve 
the statistics by binning together the new {/B-band data and 
the B-band data of B08 (lower panels of Fig. [T). The standard 
deviation is then reduced down to 1 . 1 • 10~ 5 . 

We also find a consistency between our measurements in 
the V band and those of W09, which were centered at about 
the same wavelength (Fig.[T|i. One can see that only one of the 
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Table 2 

Parameters of the RL-atmosphere. 



0.4 0.6 

Phase 

Figure 1. Stokes q and u with ±1<t error bars for HD189733b. Upper panels: 
the new UBV measurements are shown by squares and the binned B-band 
data from B08 by open circles, and the W09 measurements by crosses. The 
U and V data are shifted in vertical by ±4-l(T 4 for clarity. Lower panels: 
All the U and B data from the years 2006-2008 binned together. The mean 
error of the binned data is 1.7- 10 , and the standard deviation is 1.1-1 — 5 . 
Curves are the best-fit solutions for the Rayleigh-Lambert atmosphere. The 
normalized \ of the fit is 1.16. 

six W09 measurements (which provided the upper limit) is at 
a phase near elongation where a polarization signal of about 
5 T0~ 5 is expected. It matches perfectly the W09 upper limit. 
Therefore, the "null" result concluded by W09 is misleading 
due to the overlooked wavelength dependence of polarization. 

3. ORIGIN OF POLARIZATION 

Our primary hypothesis is that the observed polarization is 
due to scattering in the planetary atmosphere. The shape of 
the phase curve with two peaks near elongations is a strong 
argument in favor of this (Fig.[TJ. It can be reproduced with a 
simple reflecting model (Sect. |4j. 

However, one can imagine alternative scenarios for the ori- 
gin of the observed polarizati on. For examp le, starspots de- 
tected on the host star (e.g., iPont et al.ll2007l hereafter P07) 
when appear near the limb may cause incomplete cancellation 
of the stellar limb polariza tion. Assuming typ ical parameters 
of starspots on a K dwarf (Berdvugina 2005) and computing 
its mo st probable limb polarization at 400 nm ( Fluri & Stenflo 
1999), we evaluated a possible effect for the large spot Fea- 
ture A detected by P07. Such a spot near the limb would 
cause maximum 5 TO -7 polarization. If its area is artificially 
increased up to 1%, whi ch is the maximu m spot area ob- 
served on the stellar disk (Winn et al. 2007), the polarization 
can be 3T0~ 6 . Note that the projected spot area is strongly 
reduced toward the limb, and cannot be as large as 1 % near 
the limb (this is in contrast, e.g., to transits, when the pro- 
jected area of the planet does not depend on the limb dis- 
tance). Also, due to possible simultaneous presence of sev- 
eral spots near different parts of the limb, an expected effect 
is even smaller. These calculations demonstrate that the spot 
area on HD 189733 is simply too small to be responsible for 
the 10~ 4 effect in linear polarization. A transversal Zeeman 
effect due to starspot magnetic fields is also excluded at the 
level above 10" 6 dMoutou et al.ll2007l) . 

Another symmetry breaking effect resulting in a non-zero 





U 


B 


V 


RI 




10±15 


14±6 


5±20 






1.19±0.24 


1.18±0.10 


0.75±0.20 


<0.43 


PRL 


<0.62±0.30 


0.61±0.12 


0.28±0.16 


<0.09 



Note. — Q values differing by 180° are equally good. If our V band data are com- 
bined with those of W09, fi=0° ±9° , ft RL =0.79±0.13, and p=0.3 1 . 

polarization could be due to tidal interac t ion of the planet with 
the star. Following ICondon & Schmidl (119751) we have eval- 
uated the height of tidal bumps on the stellar surface to be 
about 1 km. This results in the ellipticity of the star of 2-10 -6 
and the maximum residual limb polarization of 10~ 9 , which is 
totally negligible. 

We therefore concentrate our efforts on modeling the ob- 
served polarization under the assumption that it is caused by 
scattering in the planetary atmosphere. 

4. RAYLEIGH SCATTERING IN THE PLANETARY ATMOSPHERE 

To model polarization of the light scattered off a planet we 
employ the Rayleigh-Lambert (RL) approximation, i.e. as- 
suming (i) Rayleigh scattering for polarization and (ii) the 
Lamber t sphere with the geometrical albedo pl =2/3 for in- 
tensity (Sea ger et al.l[2000t iFluri &~B erdvugina 2010, B08). 
Since the albedo of the atmosphere is fixed, the wavelength 
dependent properties of the planet are effectively included 
into the radius of the RL-atmosphere /?rl for different bands 
which represents a geometrical limit for the unity optical 
thickness in the atmosphere. Note that this model implies a 
Bond albedo of 1 and scattered light maximally polarized, so 
the values for radii and albedos should only be considered as 
limits for a spherically symmetric case. We use the same in- 
put p arameters and the y 2 min imization procedure as in B08 
and (Fluri & Berdvugina 2010) and apply it to the original 
measurements obtained in 2006-2008 for each passband sep- 
arately. For the orbit inclination i we used 94.32° as it results 
in a smaller \ 2 value as compared to the complementary value 
of 85.68° reported by P07. The best-fit values of the two free 
parameters of the model £1 and Rrl are listed in Table [2] The 
modeled Stokes q and u curves are shown in Fig.Q] The val- 
ues for the B band agree with those obtained by B08 within 
the standard deviation. If our V-band data are combined with 
those of W09, the parameters are about the same but the errors 
improve. 

To evaluate a wavelength-dependent geometrical albedo 
(Table 0, we compare Rrl with th e results by P08. They 
employed analytical transit curves by Mande l & Agoll ( 120021) 
under the assumption that the planetary body is fully opaque 
which provides lower limits at observed wavelengths. They 
modeled the transit curves by allowing the equivalent ra- 
dius of the planet, R eq , to vary with wavelength and found 
an increase of the radius by ^500 km at 0.55 /im as com- 
pared with 1 /im, and by ~ 1000 km as compared with mid- 
infrared. Using /f eq =1.1524/fj for 0.55 /urn (7?*=0.755fl Q ), 
we evaluate the geometrical V-band albedo of 0.28 by scaling 
p v = PhRl L (V)/R 2 eq (V). The transit in the B ba nd (0.44 /mi,) 
was o bserved about 0.4% deeper than in V dBouchv et alj 
12005b . which resulted in R eq = 1 .26R], but it was deduced with 
different orbital and stellar parameters than used by P08. By 
applying the P08 parameters, we obtain R eq (B) = 1 .23Rj. This 
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Figure 2. Polarimetric amplitudes Aq (top) and geometrical albedo p (bot- 
tom) measured in various passbands. Our new data are shown by open circles, 
the W09 upper limit by a cross, and the L09 upper limit for r Boo b scaled to 
the orbit of HD189733b by a triangle. Horizontal bars show the passbands' 
FWHI. Upper panel: Th in horizontal lines indicat e Aq for RjaJR] = 1.26, 
1.0, 0.8, 0.5 and f2=14° IFluri & Berdyugina 2010); solid (red) curve is the 
Rayleigh law Aq = 3.6 ■ 10 s A -4 scaled to fit the BVRI data. Lower panel: 
Solid curve is the geometrical albedo of Neptune smoothed by a 100 nm box- 
car. Downward arrows are upper limits obtained from the solid curve of the 
upper panel and the transit data by P08. Our data are consistent with other 
measurements and reveal Rayleigh scattering in the planetary atmosphere. 

leads to /?#=(). 61, which is a factor of two larger than that for 
the V-band. In the U band R eq is unknown but can be signifi- 
cantly larger tha n in the optical, as found for other hot Jupiters 
in the UV (e.g..lVidal-Madiar et alj |200l iLinskv et aTl EoTot 
iLecaveheTd es Etan gs et al.M2010lh Using R eq (B), an upper 
limit of pu is 0.62. This high albedo in the blue is consisten t 
with model predictions for hot Jupiters (Burr ows et al.l l2008). 
Using R e q(V), the upper limit for both pu and ps is 0.7. For 
red bands (RI) we can use the upper limit of the polariza- 
tion amplitude obtained by L09 for r Boo b but scaled to the 
orbit of HD189733b, which results in 1.210" 5 (Fig. |2) and 
R Rh = 0.43fij. Using R eq = 1.1469/?j for 0.85 fim from P08, 
Pri is 0.09, in agreement with uppe r limits for hot Jup iters in 
optical red wavelengths (e.g., L09, iRowe et ail 120081) . Note 
that the absolute values of R ecj and p rely on /?»=0.755/?©, 
while a direct measurement of the stellar angu lar diameter 
resulted in O.779±O.O52i? dBaines et al.ll2007l) . Using this 
value would increase radii and decrease albedos. Also, the 
geometrical albedo derived from polarimetry does not include 
unpolarized thermal emission of the planet which starts to 
dominate at longer wavelengths (e.g.. Burrows et al. 2008). 

The wavelength dependence of /?rl and p indicates 
that optical polarization is dominated by Rayleigh scatter- 
ing with a characteristic increase of the cross section to- 
ward the blue. The same conc lusion was reached by 
iLecaveher" d es Etangs et al. (2008a) who modeled the P08 
transit data with a translucent atmosphere (single layer). They 
convincingly showed that the main optical opacity source 
in the HD 189733b atmosphere must have a cross section 
with the A" 4 power law and suggested two possible scatter- 
ing sources, H2 molecules and MgSiCh condensate. The the 
radius change from the red to the blue of such an translu- 
cent atmosphere increased by a factor of two, up to about 
1000 km, as compared to the fully opaque case of P08. Note 
that extrapolation of the planet radius using this model to the 
near infrared fits well the corresponding transit measurements 
dSing et al.l 2009), but it underestimates the radius in the B 



band (and probably in U). This can be due to additional opac- 
ity sources in the optical, e.g., H" absorption which increases 
at 0.3-0.9 /im and decreases at 0.9-1.6 /xm. 

The A" 4 law fits well the BVRI polarization amplitudes 
and upper limits (Fig. |2). Thus, it appears that most op- 
tical photons (A>400nm) are scattered only once in the at- 
mosphere of HD 189733b, because in such a case the polar- 
ization am plitude is prop ortional to the scattering cross sec- 
tion (e.g., S tenflo 2005). Therefore, models with multiple 
scattering of all photons cannot explain the observed polar- 
ization in the blue as they predict the maximum polariza- 
tion for HD1897 33b of (3-4)10" 5 (e.g.. ISeager eTa!ll2000T: 
iStam et al.ll2004l L09; using the V-band transit radius). 

The U band polarization amplitude deviates from the 
Rayleigh law. Interestingly, even though the effective 
temperature of the Solar system planet s are cooler than t his 
planet, we find that p(X) for Neptune ( Karkosch kall 1 9941) is 
quite similar to HD 189733b (Fig. [2] lower panel), which is 
also similar to that of Uranus but somewhat different from 
Jupiter and Saturn. The reflected and polarized spectra of 
Neptune are well known to be strongly influenced by Rayleigh 
and Raman scatteri ng by H? as well as methane absorp - 
tion in the red (e.g., Sromovsky 20 05tlJoos & Schmidll2007l) . 
Rayleigh scattering results in the characteristic decrease of the 
albedo and polarization in the continuum at A>500nm, while 
Raman scattering reduces the reflectivity of the atmosphere at 
shorter wavelengths due to partial absorption of the energy of 
blue photons. Absorption by molecules is also possible. Thus, 
enhanced blue absorption can explain the f/-band polarization 
and albedo of HD 189733b, while Raman scattering can be an 
additional opacity source leading to heating of exoplanetary 
atmospheres. 

The atmosphere of Neptune can be approximated by a 
high-altitude haze layer above a semiinfinite cloud deck 
(Sromovsky 2005). Similarly, a prel iminary semi-em pirical 
model atmosphere of HD 189733b dBerdvuginal 1201 ll) ex- 
plains the observed polarization by the presence of a distinct 
dust condensate layer unbeneath a thin gaseous layer. The de- 
tails of this model will be presented in our forthcoming paper. 

5. CONCLUSIONS 

The new UBV polarimetric observations of HD 189733b 
confirm the first detection of polarized reflected light from the 
hot Jupiter by B08 at the 10~ 4 level with the standard deviation 
of 10~ 5 . Our data are consistent with upper limits obtained 
by others at different wavelengths. They clearly demonstrate 
the dominance of Rayleigh scattering in the planetary atmo- 
sphere at optical wavelengths longer than 400 nm. At shorter 
wavelengths an additional opacity mechanism (e.g., Raman 
scattering) plays a significant role. Our result establishes po- 
larimetry as an important tool for studying directly exoplane- 
tary atmospheres in the visible and near UV. In the near future 
it will be employed for non-transiting systems. 
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